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Full Faddeev-type calculations are performed for ^^Be breakup on proton target at 38.4, 100, 
and 200 MeV/u incident energies. The convergence of the multiple scattering expansion is inves- 
tigated. The results are compared with those of other frameworks like Distorted Wave Impulse 
Approximation that are based on an incomplete and truncated multiple scattering expansion. 
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I. INTRODUCTION 

Interest in quasi free breakup of a Halo nucleus stems 
mainly from the structure information about the wave 
function of the struck valence neutron. Qualitatively, the 
quasi free scattering is a process [l|, Q in which a particle 
knocks a nucleon out of the nucleus (A) and no further 
strong interaction occurs between the (A„_i) nucleus and 
the two outgoing particles. 

Assuming that the nucleus (A) can be well described by 
an inert (A„_i) core and a valence neutron, this reaction 
can be studied using the few-body Faddeev or equivalent 
Alt, Grassberger, and Sandhas (AGS) scattering frame- 
work [1,13, Q which can be viewed as a multiple scattering 
series in terms of the transition amplitudes for each inter- 
acting pair. Both equations are consistent with the cor- 
responding Schrodinger equation and therefore provide 
an exact description of the quantum three-body problem 
but are more suitable for the numerical solution. The 
Faddeev equations are formulated for the components of 
the wave function, while the AGS equations are integral 
equations for the transition operators that lead more di- 
rectly to the scattering observables. 

Traditionally this reaction has been analysed using the 
Distorted Wave Impulse Approximation (DWIA) scatter- 
ing framework fd\ . In this approach one assumes that the 
incoming particle collides with the struck nucleon as if it 
was free, and therefore the cross section for the quasi 
free scattering is given by the product of the transition 
amplitude between the two coUiding particles modulated 
by the wave function of the nucleon inside the nucleus. 
Furthermore the relative motion of the particles in the en- 
trance and exit channels is described by distorted waves. 
This approximation leads, in each order, to an incomplete 
and truncated multiple scattering expansion. 

In order to obtain accurate information about the wave 
function, one needs to investigate the effect of higher or- 
der rescattering between the two coliding particles and 
the (A„_i) nucleus in the calculated knockout observ- 
ables. In addition, the contribution of the neglected 
terms in each order needs to be assessed. This is the 



aim of the present work. 

We consider here the quasi free scattering of the A 
nucleus by a proton target leading to the breakup reac- 
tion p(A, A„_in)p. We take as an example the breakup 
of ^^Be halo into a proton target at 38.4, 100, and 200 
MeV/u. 

Recently the elastic scattering of ^^Be from a proton 
target was analysed using the Faddeev/ AGS multiple 
scattering approach Its was shown that even at in- 
termediate energies of 200 MeV/u one needs to take into 
account 3rd order multiple scattering contributions that 
are not included in the Glauber scattering framework Q. 

In this work we briefly review the Faddeev/ AGS mul- 
tiple scattering framework in section II. In section III 
we sumarize the DWIA formalism and identify the terms 
of the corresponding multiple scattering expansion. In 
section IV we show the calculated fivefold breakup cross 
section using the Faddeev/ AGS equations and compare 
with the results of the DWIA type. 



II. THE FADDEEV/AGS EQUATIONS 

In this section, for completeness we briefly describe 
the Faddeev/ AGS multiple scattering framework which 
treats all open channels in the same footing. Let us 
consider 3 particles (1,2,3) interacting by means of two- 
body potentials. In this section we use the odd man out 
notation appropriate for 3-body problems which means, 
for example, that the interaction between the pair (2,3) 
is denoted as vi. We assume that the system is non- 
relativistic and we write its total Hamiltonian as 



H^Ho 



7 



(1) 



with the kinetic energy operator Hq and the interaction 
for the pair 7. The Hamiltonian can be rewritten as 



where Ha is the Hamiltonian for channel a 

Ha — Ho + Va , 



(2) 



(3) 
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and represents the sum of interactions external to 
partition a 

= E "7 • (4) 

The a — partition corresponds to three free particles 
in the continuum where y° is the sum of all pair inter- 
actions. The total transition amplitude which describes 
the scattering from the initial state a to the final state (3 
is given in the post form (see Refs. [1,I3| for a derivation) 

T^" ^yl' + Y, T^Go^ (5) 

with the transition amplitude 

i-y = + VjGot.y (6) 
and the free resolvent 

Go=^iE + iO-H„)-\ (7) 

where E is total energy of the three-particle system. 
Equivalently in the prior form we write 

" = v'' + Y^ hGoTr . (8) 

One then defines the operators 11'^°' which are equivalent 
to the transition amplitudes on the energy shell, such 
that 

U^" = SfjcG-'+T^"- 

= Sp^G^^ + r!^" , (9) 

where Sfja — I - 5f3a- From Eqs. (HHH) 

= 6pa,Go' + E U^^Got^S-^c , (10) 

7 

or 

U"" = Sp^Go' + E Sp^t^GoU-^" . (11) 

7 

These are the well known three-body AGS equations Q 
which for breakup (/? = in the final state) become 

i70" = Go-^+E^7GoC/^", (12) 

7 

where 11^°" is obtained from the solution of Eq. I|lip with 
a,l3,^= (1, 2, 3). The scattering amplitudes are the ma- 
trix elements of U^°' calculated between initial and final 
states that are eigenstates of the corresponding channel 
Hamiltonian Ha [Hp) with the same energy eigenvalue. 
For elastic scattering and transfer those states are made 
up by the bound state wave function for pair a{(3) times a 
relative plane wave between particle a (/3) and pair a (/?). 
For breakup the final state is a product of two plane 
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FIG. 1: Single scattering diagrams for breakup in the Faddeev 
scattering framework. 

waves corresponding to the relative motion of three free 
particles that may be expressed in any of the relative Ja- 
cobi variables. In the latter case the contribution of the 
Gq ^ term is zero. 

The solution of the Faddeev/AGS equations can be 
found by iteration leading to 

7 7 4 




+ ••• , (13) 

where the series is summed up by the Fade method 
The successive terms of this series can be considered 
as first order (single scattering), second order (double 
scattering) and so on in the transition operators. The 
breakup series up to third order is represented diagra- 
matically in Figs. [T]-[3]where the upper particle is taken 
as particle 1 scattering from the bound state of the pair 
(23). 

In our calculations Eqs. (fTD - [T3|) are solved exactly 
after partial wave decomposition and discretization of all 
momentum variables. Since, in addition to the nuclear in- 
teraction between all three pairs, we include the Coulomb 
interaction between the proton and ^°Be, we follow the 
technical developments implemented in Refs. [l3, [lH for 
proton-deuteron and a-deuteron elastic scattering and 
breakup that were also used in Ref. [3| to study p-^'^Be 
elastic scattering. The treatment of the Coulomb interac- 
tion in the framework of three-body Faddeev/AGS equa- 
tions is based on the method of screening of the Coulomb 
interaction plus renormalization [10 . ll| ; it leads to fully 
converged results for the on-shell elastic, transfer and 
breakup observables that are independent of the choice 
of screening radius used in the calculations, as long as 
it is sufficiently large (R > 10 fm for the observables 
considered in the present work). 

III. THE DWIA/PWIA 

Let us consider the reaction A{a, ab)B where an inci- 
dent particle a knocks out a nucleon or a bound clus- 
ter b in the target nucleus A resulting in three particles 
(a, b, B) in the final state. This reaction has been ana- 
lyzed within the DWIA as described for example in the 



3 



IT 



H 



n 



IT 



n 



FIG. 2: Double scattering diagrams for breakup in the Fad- 
deev scattering framework. 
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FIG. 3: Triple scattering diagrams for breakup in the Faddeev 
scattering framework. 



work of Chant and Roos This is an approximate for- 
maHsm and its vahdity in the appUcation to radioactive 
beam studies needs to be assessed. In the work of ref- 
erence it is assumed that the projectile strucks the 
ejected particle freely and that the transition amplitude 
is given by 



Tab ^ {v^l\tab\(t'BbviA^), 



(14) 



where 0Bb is the (Bb) bound state wave function, and 
rj'^ and ry^ Jf, describe the relative motion of the particles 
in the entrance and exit channels, respectively, neglecting 
the interaction Vab- Therefore in the entrance channel 
satisfies 



(TaA + VaA " Kfc " £aA)viA = ^ 



(15) 



where CaA is the relative kinetic energy and TaA the rel- 
ative kinetic energy operator. In standard DWIA calcu- 
lations one takes VaA — Vab ~ VasiraA)- For the exit 
channel one assumes as a good approximation to write 
the 3-body wave function rf^ab ^ factorized product 



where 



and 



[Tal 



''I Bab ~ viJvbB'' ^ 



VaB - ^aBWoB 



{TbB + VbB - ebBWbV = 



0, 



(16) 



(17) 



(18) 



The potentials VaB and VbB are taken to be the optical 
potentials which describe the a + B and b + B scattering 
at energies CqB and CbB respectively. Therefore one writes 



Tab - (viJviB^\tab\Ma7)- (19) 
If we now introduce plane waves in the exit channel 

Tab {XaBXbB \{l+taBGo)il+tbBGo)tab\(l)Bbv'a^A )' (20) 

we then obtain the truncated multiple scattering series 
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tab + taBGotab + tbBGotab + taBGatbBGotab 



(21) 



The first order term is the single scattering represented 
diagramatically in Fig. |4l This is the only term retained 
in Plane Wave Impulse Approximation (PWIA) calcu- 
lations. The second and third order terms in Eq. l(2T|) 
represented in Fig. [5] are due to distortion in the exit 
channel. The relative importance of these terms to the 
scattering can be assessed by performing Faddeev multi- 
ple scattering calculations. 

The distortion in the entrance channel is more difficult 
to bridge with few-body calculations. In order to esti- 
mate this distortion in the case of a light emited particle 
one could say for instance, that 



„(+) 

VaA 



(1 + G()taB)Xa/ 



(22) 



This leads to additional multiple scattering contributions 

tabGotaB + taBG^tabGotaB + tbBGotabGotaB + 0(t'^). 

(23) 

These terms are represented in Fig. [6l We note that the 
distortion on the entrance and exit channels contribute 
to a truncated and incomplete multiple scattering series 
at each order. 

We concentrate on the study of the multiple scatter- 
ing expansion. In addition, we want to investigate if the 
single scattering term becomes the dominant contribu- 
tion as the energy increases. It was shown in the work of 
Chant and Roos Q that, contrary to what is expected, 
the DWIA results do not converge to the PWIA results 
as the projectile energy increases. The importance of 
including the distortion in comparison with the full Fad- 
deev multiple scattering series is also estimated in the 
present work. We refer as DWIA-full a multiple scat- 
tering calculation that includes the multiple scattering 
terms in Eqs. (pT]) and l(23|) . The calculations based on 
these equations follow exactly the same procedure as ex- 
plained before for the solution of the Faddeev / AGS equa- 
tions ifTTj) through lfT3|) , both in terms of initial and final 
states as well as treatment of Coulomb, so that one may 
obtain a consistent comparison between different terms 
in the multiple scattering representation of DWIA and 
Faddeev/ AGS amplitudes. 

We note that further approximations are usually made 
in standard applications of the DWIA when evaluating 
the transition amplitude, namely: 
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FIG. 4: Single scattering diagram in the DWIA scattering 
framework. 
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FIG. 5: Diagrams due to the distortion in the exit channel in 
the DWIA scattering framework. 



• The potential approximation in the entrance chan- 
nel VaA - Vab ^ VasiraA)- 

• The factorization approximation, which is only ex- 
act in PWIA. 

• On-shell approximation of the transition ampli- 
tude. 

We refer to this as DWIA-standard. An attempt to esti- 
mate the effect of such approximations was done for ex- 
ample in the work of Ref. for (p,2p) reactions. The 
Faddeev calculations and the DWIA-full do not make use 
of such approximations in the summation of the multiple 
scattering expansion, and therefore we are not address- 
ing here the validity of such approximations which are re- 
sponsible for additional shortcomings of DWIA-standard 
vis-a-vis the exact calculations. 
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FIG. 7: (Color online) Cross section for the breakup 
"Be(p,pn) at 38.4 MeV. 
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FIG. 8: (Color online) Cross section for the breakup 
"Be(p,pn) at 100 MeV. 



IV. RESULTS 

We now proceed to solve Faddeev/ AGS equations ifTD 
- [T3|) and to calculate the observables for the proton- 
^^Be breakup. As mentioned above the treatment of the 
Coulomb interaction, as well as the calculational tech- 
nique, is taken over from Refs. (lol. 

First we describe the pair interactions p-n, p-^^Be and 
n-^^Be. For p-n we take the realistic nucleon-nucleon CD 
Bonn potential • The interaction between the valence 
neutron and the ^°Be core is assumed to be L-dependent 

"n~ n~ 
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FIG. 6: Diagrams due to the distortion in the incoming chan- 
nel in the DWIA scattering framework. 



as described in [?]. For the potential between the proton 
and ^"Be core we use a phemenological optical model 
with parameters taken from the Watson global optical 
potential parametrization 0, EBl • The energy dependent 
parameters of the optical potential are taken at the pro- 
ton laboratory energy of the p-^^Be reaction in inverse 
kinematics. 

In the solution of the Faddeev equations we include n- 
p partial waves up to relative total angular momentum 
inp < 8, n-i°Be up to < 3, and p-i°Be up to L < 10. 
Three-body total angular momentum is included up to 
100. 

As a first result we show in Figs. [THS] the fivefold dif- 
ferential breakup cross section a j d£tid£t2dS versus the 
arclength S at Euh/A = 38.4 MeV, 100 MeV, and 200 
MeV respectively. The kinematical configurations are 
characterized by the polar and azimuthal angles {9i,(j)i) 
of the two detected particles labeled 1 and 2 in Fig. [101 
We assume those particles to be the ^°Be core (1) and 
the neutron (2). We take three quasi free scattering 
configurations: (conf 1) with — (0°,0°) and 
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FIG. 9: (Color online) Cross section for the breakup 
"Be(p,pn) at 200 MeV. 
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FIG. 11: (Color online) Cross section for the breakup 
"Be(p,pn) at 38.4 MeV. 
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FIG. 10: (Color online) Kinematic angles for breakup. 
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(6*2,02) = (45°, 180°); (conf 2) with (6*1, 0i) = (0°,0°) 
and (6*2,02) = (30°, 180°); (conf 3) with (6*1,01) = 
(0°,0°) and (6*2,02) = (15°, 180°). The arclength S 
is related to the energies Ei of the two detected par- 
ticles in the standard way [il as = l^dS with 
dS = [dE\ + dE"^)^!"^ . In each configuration, the soHd 
fines represent the fufi multiple scattering (MS) results 
obtained from the solution of Faddeev/AGS equations. 
The single scattering approximation (SS) represented by 
the dashed lines is calculated taking into account the pro- 
ton scattering from the struck nucleon and the scatter- 
ing from the valence core as represented by diagrams in 
Fig. [TJ As follows from Fig. [7] the SS aproximation over- 
estimates the exact result in any configuration at 38.4 
MeV/u. At 200 MeV/u, the SS already provides a good 
representation of the full calculation in the case of conf 
1, but increasingly overestimates MS in conf 2 or conf 3 
where the neutron is emitted at smaller angles. Thus, 
contrary to what is predicted by Chant and Roos [6|, if 
the neutron is not detected at very forward angles, the 
single scattering approaches the full calculation as the 
beam energy is increased. The double scattering (DS) 
approximation, represented by the circles, and calculated 
taking into account the diagrams of Figs. [1] and [21 re- 



FIG. 12: (Color online) Cross section for the breakup 
"Be(p,pn) at ICQ MeV. 



produces well the full calculation at 200 MeV/u energy. 
Therefore the Faddeev/AGS multiple scattering series for 
breakup has converged at second order level at this en- 
ergy. 

Now we compare Faddeev/AGS and DWIA-full results 
for the breakup observables. We begin by analyzing 
the single scattering term in Figs. [TlHl3l The breakup 
cross section calculated taking into account both proton- 
neutron and proton-core contributions is represented by 
the solid line. The circles correspond to the PWIA, i.e., 
without the proton-core scattering contribution. At 38.4 
MeV/u this contribution is very small in the configura- 
tions 1 and 2, but is sizable when the neutron is emmited 
in the forward region. Therefore, even at the single scat- 
tering level, calculations become innacurate in some con- 
figurations at low energies. At high energies the proton- 
core scattering contribution is negligible in all configura- 
tions. 

Since, as shown in Figs. [THSI the DS results are very 
close to the full MS results, we compare next in Figs.fTil- 
[16] the predictions of Faddeev/AGS and the DWIA-full 
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FIG. 13: (Color online) Cross section for the breakup 
"Be(p,pn) at 200 MeV. 



approach including in both calculations the correspond- 
ing multiple scattering terms up to second order. As 
explained before the DWIA-full approach provides an in- 
complete series in all orders, thought in first order the 
effects of the missing term are negleagible at high ener- 
gies as shown in Figs. [TlHTSl 

The Faddeev/AGS (solid line) results include all the 
single and double scattering terms shown in Figs 
The DWIA-full assumes that the contribution of the scat- 
tering between the proton and the neutron is dominant 
and therefore the single scattering contribution due to 
the scattering of the proton from the ^°Be core is ne- 
glected. The circles include the first and second order 
diagrams originated from the distortions in the exit chan- 
nel represented in Fig. [H and the dashed curve includes, 
in addition, the second order diagrams originated from 
the distortion in the incoming channel represented in 
Fig. [6l The DWIA-full results to second order (dashed 
line) provide a poor approximation to the Faddeev/AGS 
at 38.4 MeV in any configuration. As the energy increases 
DWIA-full approaches the exact result. However one 
should keep in mind that DWIA-standard calculations 
assume further approximations as discussed in Sect. III. 

In Figs. [9] and [TllfTSl we show that in the high en- 
ergy regime PWIA approaches the exact result at least in 
some suitable breakup configurations. Therefore, in this 
energy regime, in order to get a rehable prediction of the 
breakup observables in quasi free scattering, we advocate 
that a better approach is to work with the PWIA scat- 
tering framework within suitably chosen configurations 
rather than in DWIA where uncontrolable approxima- 
tions are made. On the contrary, in PWIA the factoriza- 
tion is exact and the transition amplitude for the scat- 
tering between the incident and knocked particle is on 
the energy shell in the weak binding limit of the knocked 
particle. 
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FIG. 14: (Color online) Cross section for the breakup 
"Be(p,pn) at 38.4 MeV. 
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FIG. 15: (Color online) Cross section for the breakup 
"Be(p,pn) at 100 MeV. 



V. CONCLUSIONS 

We have calculated the breakup cross section for ^^Be 
breakup from protons at 38.4, 100 and 200 MeV/u. Quasi 
free scattering conditions are considered in different con- 
figurations for the emitted neutron. The Faddeev/AGS 
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FIG. 16: (Color online) Cross section for the breakup 
"Be(p,pn) at 200 MeV. 
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formalism is used to show that at these energies the mul- 
tiple scattering expansion has converged at the second 
order. We have also shown that at high incident energies 
the single scattering (PWIA) approaches the full calcu- 
lation in some suitable configurations, namely when the 
neutron is not emitted at forward angles. We have also 
shown that the DWIA-fuU scattering framework provides 
a poor approximation of the Faddeev/AGS formalism at 
low energies, but approaches it at high energies. How- 
ever since standard applications make further use of un- 



controlable approximations such as the factorization and 
on-shell approximations, one advocates that at these high 
energies it is adviseable to use the exact approach or in- 
stead the PWIA in suitable configurations. 

Acknowledgements: The work of R. Crespo is 
supported by the Fundagao para a Ciencia e Tecnologia 
(FCT), through grant No. POCTI/FNU/43421/2001. 
A. Deltuva is supported by the FCT grant 
SFRH/BPD/34628/2007 and all other authors by 
the FCT grant POCTI/ISFL/2/275. 



[1] G. Jacob, TH. A.J. Maris, Rev. Mod. Phys. 38, 121 
(1966). 

[2] G. Jacob, TH. A.J. Maris, Rev. Mod. Phys. 45, 6 (1973). 

[3] L.D. Faddccv, Zh. Eksp. Theor. Fiz. 39, 1459 (1960) [Sov. 
Phys. JETP 12, 1014 (1961)]. 

[4] E.O. Alt, P. Grassberger, and W. Sandhas, Nucl. Phys. 
B 2, 167 (1967). 

[5] W. Glockle, The Quantum Mechanical Few-Body Prob- 
lem (Springer- Vcrlag, Berlin/'Hcidclbcrg, 1983). 

[6] N.S. Chant and P.G. Roos, Phys. Rev. C 15, 57 (1977). 

[7] R. Crespo, E. Cravo, A. Deltuva, M. Rodriguez-Gallardo, 
and A.C. Fonseca, Phys. Rev. C 76, 014620 (2007). 

[8] J.S. Al-Khalili, R. Crespo, R.C. Johnson, A.M. Moro, 



and I.J. Thompson, Phys. Rev. C 75, 024608 (2007). 
[9] G. A. Baker and J.L. Gammel eds.. The Fade Approxima- 
tion in Theoretical Physics (Academic, New York, 1970). 

[10] A. Deltuva, A.C. Fonseca, and P.U. Sauer, Phys. Rev. C 
71, 054005 (2005); 72, G54GG4 (2005); 73, 057001 (2005). 

[11] A. Deltuva, Phys. Rev. C 74, 064001 (2006). 

[12] N. Kanayama, Y. Kudo, H. Tsunoda, and T. Wakasugi, 
Progress of Theoretical Phys 83, 540 (1990). 

[13] R. Machlcidt, Phys. Rev. C 63, 024001 (2001). 

[14] K. Chmielewski et al., Phys. Rev. C 67, 014002 (2003). 

[15] B.A. Watson, P.P. Singh, and R.E Segel, Phys. Rev. 182, 
978 (1969). 



